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Why are accreting T Tauri stars observed to be less 
luminous in X-rays than non-accretors? 
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ABSTRACT 

Accreting T Tauri stars are observed to be less luminous in X-rays than non-accretors, 
an effect that has been detected in various star forming regions. To explain this we 
have combined, for the first time, a radiative transfer code with an accretion model 
that considers magnetic fields extrapolated from surface magnetograms obtained from 
Zeeman-Dopplcr imaging. Such fields consist of compact magnetic regions close to 
the stellar surface, with extended field lines interacting with the disc. We study the 
propagation of coronal X-rays through the magnctospherc and demonstrate that they 
are strongly absorbed by the dense gas in accretion columns. The reduction in the ob- 
served X-ray emission depends on the field geometry, which may explain why accreting 
T Tauri stars show a larger scatter in their observed X-ray luminosity compared to 
non-accreting stars. 
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1 INTRODUCTION 



Accreting T Tauri stars are observ ed to be less lumi- 
nous in X-rays than non-ac cr etors (IStelzer fc Neuhauserl 
200 ll; IFlaccomio et al.l I2003all3. [ioOrj: IStassun et all 1200? 



Preibisch et alj|2005l : iTelleschi et al.ll2007an . Accreting stars 
appear to be a factor of ~ 2 less luminous, and show a 
larger variation in their X-ray activ ity compared to non- 
accreting stars l|Preibisch et alj|2005l ). However, it is only in 
recent years that this result has become cl ear, with previ- 
ous s tud ies showing conflicting results (e.g. iFeigelson et alj 
120031 andlFl accomio et al.ll2003b| ) . The apparent discrepancy 
arose from whether stars were classified as accreting based 
on the detection of excess IR emission (a disc indicator) or 
the detection of accretion related emissio n lines. However, 
with c areful re-analysis of archival data (|Flaccomio et all 
l2003al ) and recent large X-ray s urveys like the Cha ndra 
Orion Ultradeep Project (COUP; ICetman et~ai1l2005l ) and 
the XMM-Newton Extended Survey of the Taurus Molec- 
ular Cloud (XEST; iGiidel et all |2007al ) the result is now 
clear, namely that accreting T Tauri stars are observed 
to be, on average, less luminous in X-rays than non- 
accreting stars. Although the difference is small it has been 
found consistent ly in various star forming r egions, namely 
Tauru s- Aurigae ( Stelzer &i Nc uhauscr 2001; Tcll eschi et al 



2007al), the ONC (IFlaccomio et al 



2003d ; Stassun et al 



2004 jPreibisch et al.l I2005D , N GC 2264 (IFlaccomio et al 



2003d , l2006ft and Chamaeleon I (|Flaccomio et al.ll2003ah 



It should be noted, however, that such observations 
from CCD detectors are not very sensitive to X-rays that are 
produced in accretion shocks. High resolution X-ray spec- 
troscopic measurements have indicated emission from cool 
and high density plasma, most likely associated with accre- 
tion hot spots, in several (but not all) ac creting stars (e.g. 
ITelleschi et "alll2007bl : iGunther et ai1l2007r i . In this letter we 
only consider coronal X-ray emission such as is detected by 
CCD measurements. 

It is not yet understood why accreting stars are under 
luminous in X-rays, although a few ideas have been put for- 
ward. It may be related to higher extinction due to X-ray 
absorption by circumstellar discs, h owever the COUP re - 
sults do not support this suggestion (|Preibisch et alj|2005l ). 
It may be related to magnetic braking, whereby the interac- 
tion between the magnetic field of an accreting star with its 
disc slows the stellar rotation rate leading to a weaker dy- 
namo action and therefore less X-ray emission; although the 
lack of any rotation -activity relation for T Tauri stars has 
ruled out this idea (IFlaccomio et al1l2003d : iPreibisch et al.l 
120051 : iBriggs et al] 120071 '). A third suggestion is that ac- 
cretion may alter the stellar structure affecting the mag- 
netic field generation process and therefore X-ray emission 
(jPreibisch et alj|2005l h However, the most plausible sugges- 
tion is the attenuatio n of coronal X-rays by the dense gas 
in accretion columns |Flaccomio et al. 2003c; IStassun et al.l 
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l2004lPreibisch et al.ll2005l ; lGudel et al.ll2007bf ). X-rays from 
the underlying corona may not be able to heat the material 
within accretion columns to a high enough temperature to 
emit in X-rays. Field lines which have been mass-loaded with 
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dense disc material may obscure the line-of-sight to the star 
at some rotation phases, reducing the observed X-ray emis- 
sion. In this letter we demonstrate this in a quantitative 
way by developing an accretion flow model and simulating 
the propagation of coronal X-rays through the stellar mag- 
netosphere. 



2 REALISTIC MAGNETIC FIELDS 

In order to model the coronae of T Tauri stars we need 
to assume something about the form of the magnetic field. 
Observations suggest it is compact and inhomogeneous and 
may vary not only with time on each star, but also from 
one star to the next. To capture this behaviour, we use 
as examples the field structures of two different main se- 
quence stars, LQ Hya and AB Dor det ermined from Zeeman- 
Doppler imaging (|Donati et al.ll2003h . Although we cannot 
be certain whether or not the magnetic field structures ex- 
trapolated from surface magnetograms of young main se- 
quence stars do represent the magnetically confined coro- 
nae of T Tauri stars, they do satisfy the currently avail- 
able observational constraints. In future it will be possible 
to use real T Tauri magnetograms derived from Zeeman- 
Doppler images obtained using the ES PaDOnS instrumen t 
at the Canada- France-Hawaii telescope (|Donati et al.|[2007h . 
However, in the meantime, the example field geometries 
used in this letter (see Fig. [TJ capture the essential fea- 
tures of T Tauri coronae. They reproduce X-ray emission 
measures (EMs) a nd coronal densitie s which are typical 
of T Tauri stars (| Jardine et al.l |2006h ■ The surface field 
structures are complex, consistent wit h polarisation mea- 
surements (jValenti fc Johns-Krullll2004 ) and X-ray emitting 
plasma is confined within unevenly distributed magnetic 
structures close to the stellar surface, giving rise to signifi- 
cant ro tational modulation of X-ray emission (|Gregorv et al.l 
l2006bh . 

We extrapolate from surface magnetograms by assum- 
ing that the magnetic field B is potential such that V x B — 
0. This condition is satisfied by writing the field in terms of 
a scalar flux function <]/, such that B = — V\P. Thus, in or- 
der to ensure that the field is divergence- free (V.S = 0), 9 
must satisfy Laplace's equation, V 2 '!' = 0; the solution of 
which is a linear combination of spherical harmonics, 



i m r 



(1) 



where Pi m denote the associated Legendre functions. The 
coefficients ai m and bi m are determined from the radial 
field at the stellar surface obtained from Zeeman-Doppler 
maps and also by assuming that at some height R s above 
the surface (known as the source surface) the field be- 
comes radial and hence Bg(R 3 ) = 0, emulating the effect 
of the corona blowing open fi eld lines to form a stellar wind 
|Altschuler fc Newkirkl [l96 9). In order to extrapolate the 
field we used a modified version of a code originall y devel- 
oped bv lvan Ballegooiien. Cartledge fc Priestl (|l99S 



2.1 The coronal field 

For a given surface magnetogram we calculate the extent 
of the closed corona for a specified set of stell ar parameters. 
As thi s pro cess has been d e scribed in detail by I Jardine et al.l 
|2006r i and lGregorv et all (|2006al lbh we provide only a brief 
outline here. We assume that the corona is isothermal and 
that plasma along field line loops is in hydrostatic equilib- 
rium. The pressure is calculated along the path of field lines 
loops and is set to zero for open field lines and for field 
lines where, at some point along the loop, the gas pres- 
sure exceeds the magnetic pressure. The pressure along a 
field line scales with the pressure at its foot point, and we 
assume that this scales with the magnetic pressure (po = 
KB 2 ). This technique has been used successfully to cal- 
culate mean coronal densities and X-ray EMs for the Sun 
and other main sequence stars (I Jardine et al.ll2002l ) as well 



as T Tauri stars (| Jardine et al.l 120061 ). The extent of the 
corona depends both on the value of K and also on Bo 
which is determined directly from s urface magnet o grams 
and varies across the stellar surface. Ijardine et al.l (|2006h 
determined the best value of K for a given surface mag- 
netogram which results in the best fit to the X-ray EMs 
of stars from the COUP database. We have adopted the 
same values in this letter. We note that we make a con- 
servative estimate of the location of the source surface by 
calculating the largest radial distance at which a dipole 
field line would remain closed, with the same average field 
strength. The AB Dor-like coronal field has an X-ray ElvQ of 
log EM = 53.73 cm -3 (without considering accretion) and 
a mean EM-weighted coronal density of logn = 10.57 cm -3 , 
consis tent with estimate s from the modelling of individual 
flares jFavata et al.ll2005T ). The LQ Hya-like field has a more 
extended corona and consequently a lower coronal density 
and EM, log EM = 52.61 cm -3 , log n = 9.79 cm" 3 . 



2.2 The accreting field 

We consider a thin disc and assume that the disc normal 
is aligned with the stellar rotation axis. We consider a 
steady state flow model and assume that the structure of 
the magnetic field remains undistorted by the in-falling ma- 
terial and tha t the m agnet osphere rotate s as a so lid body. 
Ijardine et all (|2006t ) and iGregorv et""aH l|2006ah demon- 
strated that for most stars the closed field region is confined 
close to the stellar surface. T Tauri stars have hot coronae, 
with temperatures in excess o f several tens o f mega -Kelvin 
[e.g. see the COUP database. Ibetman et all (|2005t )]. Such 
high temperatures lead to closed field lines being opened 
by the gas pressure, limiting the closed field region to close 
to the stellar surface. In s uch cases the inner di sc may sit 
in a reservoir of open field (|Gregorv et al.|[2006aT ), This can 
be seen in Fig. [T] where the field lines that are carrying 
accreting gas are open. The question of where the disc is 
truncated remains a major problem for accretion models. It 
is still unknown if the disc is truncated in the vicinity of the 
corotation radius, the assumption of traditional accretion 



1 The X-ray EM is given by EM = J n 2 dV where n and V 
are the coronal density and volume. The EM-weighted density is 
n = J n 3 dV/ J n 2 dV. 
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Figure 1. The model coronal (first row) T Tauri magnetic fields extrapolated from the AB Dor (left-hand column) and LQ Hya (right- 
hand column) surface magnetograms, with the corresponding X-ray corona (second row) assuming a stellar inclination of 60°. Also shown 
is the structure of the accreting field (third row) and the X-ray emission images assuming that accretion is taking place (fourth row) - 
notice the occulted X-ray bright regions. For some lines-of-sight the X-ray bright regions are lightly obscured by the accretion columns, 
reducing the observed X-ray emission. For other lines-of-sight the coronal X-rays cannot penetrate the dense accreting gas. The average 
reduction in the observed X-ray EM across an entire rotation cycle is factor of 1.4 (2.0) for the AB Dor-like (LQ Hya-likc) field. The 
images are not to scale. 



models (e.g. iKonigll 11991b . or whether it extends closer 
to the stellar surface (e.g. iMatt fc Pudritj|2005l ). For this 
letter we consider that accretion occurs over a range of 
radii within the corotation ra dius. This is equivalent to the 
appro ach tak en previously by Hartmann. Hewett fc Calvetl 

Jl994l). iMuzeroile. Calvet fc Hartmannl ll200ll) 

ISvmington. Harries fc Kurosawa! (|2005l ) and [Azevedo et al.l 
|2006l ) who have demonstrated that such an assumption 
reproduces observed spectral line profiles and variability. 
It should also be noted that the accreting field geometries 
which we consider here are only snap-shots in time, and 
in reality will evolve due to the interaction with the disc. 
The accretion filling factors are 0.9% for the AB Dor-like 
field and 1.2% for the LQ Hya-like field, smaller than would 



be expected for accretion to a dipo le, but consistent with 
obser vationally inferred values (e.g. IValenti fc Johns- Krull 
12004 ). 

We assume that material is supplied by the disc and 
accretes onto the star at a constant rate. For a dipolar mag- 
netic field accretion occurs into two rings in opposite hemi- 
spheres centred on the poles. In this case, half of the mass 
supplied by the disc accretes into each hemisphere. For com- 
plex magnetic fields accretion occurs into discrete hotspot s 
distributed in latitude and longitude (|Gregorv et alj|2006af ) . 
It is therefore not clear how much of the available mass from 
the disc accretes into each hot spot. We use a spherical grid 
and assume that each grid cell within the disc which is ac- 
creting supplies a mass accretion rate that is proportional 
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to its surface area. For example, if an accreting grid cell has 
a surface area that is 1% of the total area of all accreting 
grid cells, then this grid cell is assumed to carry 1% of the 
total mass that is supplied by the disc. Therefore, as an ex- 
ample, if grid cells which constitute half of the total area 
of all accreting cells in the disc carry material into a single 
hotspot, then half of the mass accretion rate is carried from 
the disc to this hot spot. In this way the accretion rate into 
each hot spot is different and depends on the structure of 
the magnetic field connecting the star to the disc. 



2.3 Accretion flow model 

We consider a star of mass 0.5 Mq, radius 2R©, rotation 
period 6 d, a coronal temperature of 20 MK and assume that 
the disc supplies a mass accretion rate of 10 -7 Mq yr _1 . In 
order to model the propagation of coronal X-rays through 
the magnetosphere we first need to determine the density 
of gas within accretion columns. Provided that the flow of 
accreting material is along the path of the field (v is parallel 
to B) then 

^ = const. (2) 
B 



where p is the mass density (see e.g. iMestell [l968h ■ If an 
individual accreting field line carries material onto an area 
of the stellar surface of A* with velocity «* and density p*, 
then the mass accretion rate into A* is, 



M = p*v*A*. 



(3) 



It then follows from ^ and © that the density profile along 
the path of the field line may be written as, 



p(r) 



B{r) M 
B(R*) A*v(r)' 



(4) 



The density profiles therefore do not depend on the absolute 
field strength, but instead on how the field strength varies 
with height above the star. The density profiles are typically 
steeper than those derived for accretion flows along dipolar 
field lines since the strength of a higher order field drops 
faster with height above the star (see Fig. [2]). Gas is assumed 
to free-fall under gravity along the path of the field after 
leaving the disc with a small initial velocity of lOkms 
|Muzerolle et al.l 120011 ). Accretion occurs from a range of 
radii as discussed in §2.2. For the AB Dor-like field accretion 
occurs from ~ 2.0 — 2.5i?», whilst for the LQ Hya-like field, 
which is more extended, from ~ 2.4 — 2.97?*. 

In order to estimate the number density of gas within 
the accretion columns we require the temperature distri- 
bution along the path of the field lines. The temperature 
variation within the accretion columns of T Tauri stars is 
uncert ain. The most comprehensive model is that of lMartinl 
(|l996t ) who considers both heating and cooling processes. 
Here we elect to take the simplest approach and assume 
an isothermal accretion flow temperature of 7000 K - a rea- 
sonable estimate for the temperature of ac creting material 
(see e.g. lHartmann et al.|[l994l : lMartinlll996l ) ■ We obtain the 
number density profiles from, 



n(r) = p(r)/p,m H , 



(5) 



G 
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Figure 2. Some example density profiles (solid lines) for ac- 
cretion along a small selection of the complex field lines shown 
in Fig. [T] (lower left panel) assuming a mass accretion rate of 
10 — 7 Mq yr _1 . Also shown for comparison are the density pro- 
files for accretion along dipolar field lines (dashed lines), r is the 
spherical radius. 



the number density along the paths of a selection of ac- 
creting field lines, with those obtained for dipolar field lines 
shown for comparison. For our assumed accretion rate of 
lO^Mgyr -1 the flow densities range from log?i ~ 12 — 
14 cm' 3 , whilst for a lower accretion rate of 10 -8 Mq yr" 1 



where run is the mass of a hydrogen atom and p the di- 
mensionless atomic weight. Fig. [2] shows the variation of 



the range is logn ~ 11 — 13 cm 



3 SIMULATED X-RAY VARIABILITY 

We extend our previous work on optically thi n X-ray emis- 
sion from t he coronae of young stars (e.g. I Jardine et al.l 
120021 . l200fj| ; I Gregory et ail l2006bf ) to include the effects of 
absorption by the dense accretion columns. In this model 
we assume the X-ray emission from the 20 MK coronae is 
optically thin, but that the X-rays may be subsequently ab- 
sorbed in the cool (7000 K) and hence optically thick accre- 
tion columns. For the X-ray absorptive opacity we adopt a 
value of a — 10~ 22 cm 2 H _1 , typical of neutral gas at tem- 
peratures below 10 4 K a t X-ray energies of a few keV (e.g. 
iKrolik fc Kallmanlll98i ). At these energies the opacity of 
hot gas (above 10 7 K) is sev eral orders of magnitude lower 
(e.g. lKrolik fc Kallmanlll984l , Fig. 1) justifying our assump- 
tion that the 20 MK coronal X-ray emission is optically thin. 

As in our previous work, for the radiation transfer we 
use Monte Carlo techniques and discreti se the emissivity and 
density onto a spher ical polar grid (e.g. Ijardine et al]|2002l ; 
IWhitnev et al.l 120031 ) . In our simulations we adopt a stel- 
lar inclination of i = 60°. For inclinations above i > 75° the 
star and hence coronal emission, will be blocked by the verti- 
cally extended (i.e. flared) discs typical of T Tauri stars (e.g. 
iD'Alessio et~al]|l998l ; IWood et al.ll2002l ). Because we are fo- 
cusing on inclinations where the star is not blocked by the 
disc, our simulations only require our emissivity and density 
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grid to include the accretion columns and therefore our grid 
extends to the inner edge of the disc. In the Monte Carlo X- 
ray radiation transfer simulations we assume the scattering 
opacity is negligible, so our results in Fig. [T] show the ef- 
fects of attenuation of the coronal emission by the accretion 
columns. The second row in Fig. [T] shows the X-ray images 
in the absence of attenuation (i.e. X-ray opacity in the ac- 
cretion columns is set to zero) whilst the fourth row shows 
the same X-ray emission models, but with our adopted value 
for the soft X-ray opacity in the accretion columns. 

For the AB Dor-like field structure the observed X-ray 
EM is reduced by a factor of 1.4 when accretion flows are 
considered, whereas for the LQ Hya-like field the reduction 
is by a factor of 2.0 (where the reduction factor is the aver- 
age for an entire rotation cycle). For the AB Dor-like field 
there are large accretion curtains which cross the observers 
line-of-sight to the star as it rotates (see Fig.rjJ). For the LQ 
Hya-like field accretion is predominantly along field lines 
which carry material into low latitude hot spots, however, 
one of the brightest X-ray emitting regions is obscured by an 
accretion column which attenuates the coronal X-rays and 
produces a large reduction in the observed X-ray emission. 
This immediately suggests that the geometry of the accret- 
ing field is a contributory factor in causing the large scatter 
seen in the X-ray luminosities of accreting stars. 



4 CONCLUSIONS 

We have demonstrated that the suppression of X-ray emis- 
sion in accreting stars apparent from CCD observations can, 
at least in part, be explained by the attenuation of coro- 
nal X-rays by the dense material in accretion columns. The 
reduction in the observed X-ray emission depends on the 
structure of the accreting field. For stars where accretion 
columns rotate across the line-of-sight, X-rays from the un- 
derlying corona are strongly absorbed by the accreting gas 
which reduces the observed X-ray emission. This however 
does not rule out the fact that other mechanisms may also 
be respo nsible for reducing the X-ray emission in accret- 
ing stars, [jardine et al.l (|2006h have demonstrated that some 
stars (typically those of lower mass) have their outer coronae 
stripped away via the interaction with a disc. This also re- 
duces the observed X-ray emission and this effect, combined 
with the radiative transfer calculations presented in this let- 
ter, is likely to lead to a larger reduction in the observed 
X-ray emission. This would reduce the number of field lines 
which co uld be filled with coro nal gas, such as is also sug - 
gested bv lPreibisch et al.l (|2005h and lTelleschi et all (|2007ah . 
with the observed X-ray emission being further reduced due 
to obscuration by the accreting gas. A subsequent paper will 
give a detailed analysis of both of these effects drawing on 
X-ray data from the COUP and XEST. 
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